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ABSTRACT 

We present an analytical model for light echoes (LEs) coming from circumstellar material (CSM) 
around Type la Supernovae (SNe la). Using this model we find two spectral signatures at 4100 A 
and 6200 A that are useful to identify LEs during the Lira law phase (between 35 and 80 days after 
maximum light) coming from nearby CSM at distances of 0.01-0.25 pc. We analyze a sample of 89 SNe 
la divided in two groups according to their B — V decline rate during the Lira law phase, and search 
for LEs from CSM interaction in the group of SNe with steeper slopes by comparing their spectra with 
our LE model. We find that a model with LEs -I- pure extinction from interstellar material (ISM) fits 
better the observed spectra than a pure ISM extinction model that is constant in time, but we find 
that a decreasing extinction alone explains better the observations without the need of LEs, possibly 
implying dust sublimation due to the radiation from the SN. 


1. INTRODUCTION 

Type la supernovae (SNe la) are one of the most stud¬ 
ied objects in astronomy. Owing to their standardizable 
high luminosities that make them un rivaled distance in - 
dicators up to high redshifts {z ^ 2, | Jones et al.||2013|), 
astronomers searching for the ultimate fate of the uni¬ 
verse have studied them in ever greater detail, di scov¬ 


ering sever al thousands of them up to now (e. g. Sako 
et al. [2014 1. Despite the success as cosmological probes, 
increasing observations and clues reveal that the puzzle 
of the nature and mechanism generating these colossal 
explosions is still far from reaching a conclusive solution. 

As best shown in recent well-studied nearby objects 
(Nugent et al. 2011 |Bloom et al. 20121, SNe la may 
originate from the explosion of a compact CO-rich white 
dwarf (WD) in a binary system. A common candidate 
considered for the binary companion has for long been 
a non-degenerate star such as a main-sequence or a red 
giant star, that donates mass to the WD (single degen¬ 
erate, SD, scenario) either in a stable fashion so that 
the WD nears the Chandrasekhar mass and explodes 
first subsonically and then superson i cally (Ch-SD, e. g, . 
Blondin et al. 2013 Sim et al. 2013 Ropke et al. 20121 
or via unstable accretion leading to ah initial detonation 
in the outer layer of the WD that triggers a subsequent 
detonation near the center prio r to reaching the Chan¬ 
drasekhar m ass (sub Ch-SD, e. g. Sim et al. 2012 Kromer 

eTalllMol ). 


From an observational point of view, the SD scenario 
model has both evidence for and against. Among the 
observations that disfavor the model are: the absence of 


hydrogen and helium in their spectra (Lundqvist et al. 
2013[|Shappee et al.| 2013|); as well as the absence of radio 
and A-ray emission ([Horesh et al. 2012 Chomiuk et al. 


20121 which set tight constraints on mass loss from a pro¬ 


genitor; the non-detection of early emission from shock 
interaction with a companion (Bloom et al.|2012 Biancoj 


et ^ 2011 Ha yden et a l. 


non-detections (Li et al. 
red giants and i 


20 


WT 


0|) and the pre-explosion 

that generally rule out 
e-stars; tJie lack of sufficient galactic X- 


ray emission (iGilfa nov fc Bogdan|2010[ [Di Stefano|2010| 
and UV radiation ((Woods & Cntanov||2013t (Johansson 


et al.|[2M4) expected from mass accretion m the Ch-SD 


scenario; and the measured SN la rate as a function of 
redshift which challenges the modelled de lay time distri¬ 


bution for the classical Ch-SD scenario (Maoz & Man- 
i]2012 and references therein). Some of these issues 


nucci 


can be addressed through different alter native SD models 
( Justham||2011 Di Stefano et al.||2011). 

Among the observational evidence that favors the SD 
channel is the presence of nearby circumstellar material 
(CSM), presumably from mass loss in the progenitor sys¬ 
tem prior to explosion, manifested from: CSM/ejecta in - 


teraction |Hamuy et al. 2003 


and m ost notably in the case o: 


Silverman et al. 2013c) 


F'l'l'llkx (Dilday et al. 


20121; the discovery of narrow Na 1 D absorption lines 


that vary with time I Pata^^et^ []200H Pl ^din et al. 


2009 Simon et al.||20d9[ Sternberg et al.||2013p ; and the 


statistical preference tor the interstellar lines to show 


blueshifts (Sternber g et al. 2011b Maguire et al. 2013 
Phillips et al.|2013r Altogether, the observational ev¬ 


idence suggests the possibility of multiple channels for 
SNe la. 

Additionally, it has been suggested that such nearby 
CSM could affect the colors of SNe la through light scat¬ 
tering in the line of sight and explain in this way some 
of the differences in total to selective extinction ratios 
(Rv) found in SN hosts compa. r ed to the Milky Way, 
MW (Wang 2005 Goobar 2008 Amanullah & Goobar 
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20111. Heavily extincted SNe clearly show a different Ry 
while SNe with moderate extinction show values consis- 


tent with the MW (e.g. 

Mandel et al.||2011 

Burns et al. 

2U14I. Cosmological studies using sta 
fitters obtain a luminosity-color relatioi 
dening laws lower than the MW (e.g 
however it is possible that SN intrins 

mdart 
r that 
Guy 

C COl( 

iigh1 

et al. 

Drs ar 

:-curve 
st red- 
2005), 
3 more 


complicated and incorrectlv modeled (e.g. Scolnic et al. 

2014 

Chotard et al.||2011 

Conley et al.||2007). 


colors not only affords the opportunity to understand 
their nature, but also to remove a major source of sys- 
tematic uncertai nty in SN la cosmology. As shown in 
(20131, color might indeed be the principal pa- 
rameter of diversity in SN la light-curves, followed only 
by the well-known lig ht-curve width parameter of SNe la. 


Forster et al. (2013 hereafter F13) showed that the SN 
la color evolution with time is related to the strength of 
the narrow sodium absorption, suggesting that at least 
some part of it might originate from a closer interaction 
with dust than with host interstellar medium (ISM). In 
particular, they found that redder objects at maximum 
light have stronger narrow absorption lines and evolve 
faster from red to blue during the late time evolution of 
the Lira law decline of 30 to 90 d ays past maximum light 
(Lira 1995 |Phillips et al. 19991. Possible explanations 
tor this are light-echoes from (JSM that affect the late¬ 
time colors or, alternatively, CSM dust sublimation in 
the line of sight. 

In this paper, we aim to investigate the results of F13 
further and test the hypothesis of nearby CSM by look¬ 
ing for spectroscopic signatures of light-echoes in a large 
sample of nearby SN la spectra in the Lira law phase. 
Echoes from nearby ISM (and possible CSM) have previ- 


SN1992G and SN2014J at nebular phases dSchmidt et al. 

19941 ICappellaro et al. 1120011 

Wang et al.p2008a[ (Silver- 

man et a .j2013a Crotts||2014 

). These Individual studies 


SN, tens to hundreds of parsecs away, scattering hun¬ 
dreds of days past maximum light. 

Here we search for light echoes (LEs) at earlier times 
(> 30d past maximum light), coming from nearby CSM 
dust that is at less than a parsec from the SN. Such CSM 
can p otentially affect the colors an d the Lira law decline 
rate (Amanullah & Goobar 2011|). Hence we search for 
LEs in the group ol SJNe analyzed in F13 that presented 
more extinction and a steeper than normal B — V evolu¬ 
tion (hereafter fast Lira decliners) to test the hypothesis 
that these may originate in regions of nearby CSM. To do 
this, we use SNe with a shallow Lira law slope (hereafter 
slow Lira decliners) as a reference set of SNe without 
CSM interaction. 

In l|^we present our simple LE model. Then in (Owe 
focus on the prediction of observable spectroscopj^ea- 
tures to look for li^t echoes. In (|^we present the data 
and the analysis, (^summarizes the results of our search 
for LEs and in (|^we discuss the success of our LE model, 
the validity of our assumptions and other possible mech¬ 
anisms that could explain fast Lira decliners. Finally, 
the main conclusions are summarized in 0 


2. LIGHT ECHO MODEL 


The effects of the interaction of light with intervening 
dust from CSM causing scattering away from the line of 
sight, and therefore extinction and reddening, has been 
studied and modeled in depth in the past (e. g. Chevalieii 


1986 |Wang|[2005| |Patat|2005[ [Patat et al.||2(JL)D[ |Goobar| 


20081 ! We present here a simple analytical model that 


IS easy to implement numerically and makes clear ob¬ 
servable predictions to directly compare with data. The 
CSM consists of a simple spherically isotropic shell of 
dust (with Ry = 3.1), that absorbs and scatters the light 
of the SN. The radius of the shell is initially fixed at 0.05 
pc to produce LEs reaching the observer with time delays 
of ~50 days and affecting the colors during the Lira law 
phase, as light emitted at maximum light is observed at 
later epochs. At these distances we expect the tempera¬ 
ture of the dust to be slightly lower than the sublimation 
temperature (2,000 K). We only consider single scatter¬ 
ing for simplicity and because multiple scattering in the 
CSM becom es importan t when its optical depth is larger 
than 1 (see Patat| 2005|). According to our analysis all 
the SNe we considered have a total optical depth < 1 in 
the visible, with the exception of SN1997cw, SN1999gd, 
SN2003cg and SN2006X which have a total extinction 
Ay > 1 (see (j4.3|); but their Ay due to the CSM ex- 
tinct ion is lower than unity according to our models (see 
(5.2). Hence, ignoring multiple scattering is a reasonable 


approximation. 

We assume that the observed flux is the sum of the light 
coming directly from the SN and the SN light scattered 
by the CSM, i.e. light echoes. The direct flux contri¬ 
bution from the SN at a given epoch t that is extincted 
and scattered by intervening dust without including the 
contribution from LEs, can be written as: 


fit, A) = fit, (1) 

where f is the intrinsic flux of the SN, (Ts(A) and tTa(A) 
are the cross sections of the dust particles for scattering 
and absorption of photons at wavelength A, respectively. 
Ntot is the total column density of dust between the SN 
and the observer and is equal to A^ism + Acsm- Then, 
the total observed flux adding the LE contribution can 
be expressed as: 

F(CA) = /(t,A)+LE 

= fit, -f Sit, ( 2 ) 

where the first term is the intrinsic flux of the SN ex¬ 
tincted by the total column density of dust (eq. and 
the second term is the contribution from LEs, Sqt, A), 
extincted by ISM dust. cr(A) is the sum of the scattering 
and absorption cross section. We can express S'(t, A) as 
the sum of the light scattered by the CSM at different 
angles and epochs: 

I /“(ir,A)tu(A)(l-e-'™"™)*(9,A)dn (3) 

(4) 

c c 

where w(A)(l — 6“'^^^^'^'^®“) represents the fraction of 
light scattered by the CSM, which has a column density 
of A^csm- c is the speed of light, w(A) is the dust albedo 
and tr is a pseudo retarded time, i.e. the time at which 
a scattered pulse of light has to be emitted to reach the 
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observer at the same time than a pulse emitted at time t 
going straight to the observer. D is the distance between 
the photosphere of the SN and the observer and D'{9) is 
the path length travelled by a photon being scattered by 
the CSM to the observer at an angle 9 (see figure [^. 





Fig. 1. — Sketch of the CSM model. Part of the light emitted 
by the SN is scattered by the CSM and redirected to the observer 
arriving as a light echo delay. 


$(6*, A) is the scattering phase function proposed in 


Henyey & Greenstein (19411: 


<^{ 9 ,\) 


1 i-g(A)^ 

47r (1 + g{\Y — 2g{\)cos[9)Y/'^ 


(5) 


where g is the degree of forward scattering. When g = 1 
we have complete forward scattering and g = 0 means 
isotropic scattering ($(0,A) = l/(47r)). We also define 
the parameter /csm = VcsM/-Vtot to simplify the no¬ 
tation, taking values between 0 and 1. Finally, dehn- 
ing a delay parameter t = t — tr, an extinction factor 
V(A) = and making a change of vari¬ 

able, using equation]^ we obtain 

^ A)dT (6) 

where r^ax = 2i?/c is the maximum delay for a light echo 
(9 = tt) and 4’'(r, A) = 47r$(0(r), A). We performed this 
integral numerically using the Simpson’s l/3rd integra¬ 
tion rule and a time step of one day to simulate spectra 
and to fit this model to real data. 


3. LIGHT ECHO MODEL PREDICTIONS 

To simulate spectra with different extinctions and LEs 
using eq. [^we need to adopt a dust albedo, an extinc¬ 
tion law, a phase function or g(A), and spectral templates 
with no extinction nor LEs at different epochs. For the 
extinction law we take the parametrization proposed in 
Fitzpatrick (1999). We use the albedo u>(A) and the de- 
gree of forwar d scattering g(A) from the MW used in 
Goobar ([2008 ) which accounts for the dust properties of 
the GSM. We construct unreddened spectral templates 
at different epochs from weighted bootstrapped averages 
of obse rved spectra of slow B — V Lira law decliners (see 
section 4.2), together with light-curve templates that we 


need since we no rmalize the ^ectra by their V-band flux 
(see section 4.1). In Figure!^ we show the different sce¬ 
narios for late-time (Lira law phase) model spectra when 
pure extinction and simulated LEs affect the SN emis¬ 


sion. 


We search for a way to distinguish if part of the dust 
found at maximum light is producing light echoes. LE 


spectra are integrated spectra weighted by the light- 
curve, and thus dominated by spectra around peak (see 
Figure with peak template spectrum and LE spec¬ 
trum). LE spectra are blue and have very strong broad 
emission and absorption lines, with prominent peaks at 
4000, 4600, 4900 A and minima at 4400 and 6200 A. 

When LE spectra are added to SN spectra: (1) the fact 
that the LE spectra are blue has a low-order effect on 
the observed spectra by making the colors bluer, similar 
to less reddening and thus difficult to differentiate; (2) 
the strong broad lines add an additional modulation to 
the observed spectra that is very distinct to the effect 
of reddening, since it introduces differences on scales of 
a couple of hundred armstrongs. By looking specifically 
at the wavelengths where the LE spectra has peaks or 
minima, it is possible to differentiate between the two 
scenarios. In § [^ we compare these simulations to the 
observed spectra of SNe la. 

In this simulation the main signature due to LEs is 
found near 4100 A. This can be seen in Figure where 
the shape of the spectrum gets considerably mtmified in 
the LE scenario (purple and blue lines) producing a char¬ 
acteristic signature. On the other hand, in the pure ex¬ 
tinction scenario if the column density is reduced (black 
line), it produces just a smooth change in the spectrum 
compared to the same spectrum with extinction (reddest 
line). In particular, the shape of the feature near 4100 A 
will not be affected. 




Fig. 2. — Upper panel: simulated spectra at 50 days past maxi¬ 
mum light extincted by different amounts of dust (A\a=0-0.5) with 
the same Ry=3.1 extinction law. Ah spectra have been normal¬ 
ized to the same V-band flux. Bottom panel: simulated spectra 
at 50 days past maximum light with the same reddening law, but 
also with LEs due to CSM are shown. We fixed Ay = 0.5 and 
varied the fractional amount of CSM (/csm) with a radius of 0.05 
pc. The green line represents a typical maximum light spectrum, 
while in orange the LE spectrum is shown on an arbitrary scale. 
All models are normalized to the same U-band flux. 


These features change with the distance R between 
the CSM and the SN. Reducing the distance is analo- 












































4 



Fig. 3.— Simulated spectra at 50 days past maximum light with 
total extinction Ay = 0.5 and liy = 3.1, and LEs due to CSM. 
The different lines represent different fractional amounts of CSM 
(/csm) with a radius of 0.05 pc and a fixed total column density. 
The black line represents the original spectrum at 50 days past 
maximum light without extinction or LEs, while in orange the LE 
spectrum is shown in an arbitrary scale. All models are normalized 
to V-band flux. 



Fig. 4.— Features considered in the analysis as possible indicators 
of LE signatures: Ca II H&K complex at 3500-4100 A, “line 1” at 
3800-4400 A, “line 2” at 4250-4800 A and “line 3” at 3800-4800 A, 
the union of line 1 and line 2. Lines 2 and 3 come from a blend of 
Mg II and Fe II. The shaded regions indicate the pW for Ca IIH&;K, 
lines 1 and 2 and the dotted line shows the pseudo-continuum of 
line 3. 


gous to making tr = t (see eq. obtaining the result 
of a pure extinction model where S(t,X) = w(A)(l — 
On the other hand, increasing R makes 
the light echoes more diluted in time reducing the ratio 
between the LE and the intrinsic flux /. Hence, there is 
an optimum R near 0.05 pc in which LEs can affect the 
colors and spectra during the Lira phase. 

To further investigate the LE effects on spectra, we fo¬ 
cus on the range between 3000-5000 A, particularly on 
the absorption lines and their related measurable quan¬ 
tities such as the equivalent width and the slope of the 
continuum. We define five characteristic features, four 
of them that are presented in Figure]^ the Ca II H&K 
complex between 3500-4100 A, and two absorption fea¬ 
tures originating mainly from a blend of Mg II and Fe II 
at 3800-4400 A (“line 1”) and 4250-4800 A (“line 2”); 


these two together form another larger feature at 3 800- 
4800 A (“line 3”) , equivalent to the feature pW3 in 


Fo- 


latelli et al. (2013). Additionally, we use the line feature 


around the Si II absorption around 5800-6300A (“line 
4”) 

For all LE and extinction models we measure pseudo¬ 
equivalent widths ( pW), i.e. with a pseud o -continuum in 


a simi l ar fashion to Garavini e t al. (2007b); Bronder et al. 

|(|2013|). We use a semi-automatic 
IS for the 


(2008); Folatelli etTuT 


algorithm that searches for the pseudo-continuum at spe¬ 
cific regions defined in Table [T The algorithm also cal¬ 
culates the slope of the pseudo-continuum, which we find 
to be another good indicator of LEs. In figure [5 we show 
the predicted evolution of the pWs of line 1 for a SN 
with pure extinction (green lines) and with LEs due to 
CSM (black lines). The pWs vary significantly during the 
Lira phase for LE models with different CSM fractions, 
/csm, whereas they do not for different amounts of ex¬ 
tinction. This effect is strongest at larger CSM distances 
of 0.05-0.25pc. The difference between the scenario with 
extinction and with LEs is still clearer and less unbiased 
by the SN intrinsic pWs if we normalize by the pW at 
maximum light (hereafter pW ratio). 


From the five pW ratios and respective slopes, we find 
that the best candidate to be a CSM indicator is line 
1, the other lines show less differences in their evolution 
between different extinctions and CSM scenarios. 

TABLE 1 

Feature definitions for pseudo-equivalent width and 

PSEUDO-CONTINUUM. 


Feature 

Blueward limit range (A) 

Redward limit range (A) 

Ca II 

3500 - 3800 

3900 - 4100 

line 1 

3800 - 4100 

4250 - 4400 

line 2 

4250 - 4400 

4400 - 4800 

line 3 

3800 - 4100 

4400 - 4800 

line 4 

5800 - 6000 

6100 - 6300 



Days after maximum 


Fig. 5.— Pseudo-equivalent width evolution of the line 1 feature 
for simulated spectra. The green lines represent the evolution of 
line 1 with different amounts of extinction (Ay), while in black 
with different CSM amounts (/csm)? with flxed radius of 0.05 pc 
and a total column density or Ay = 1.0. 
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4. COMPARISON WITH OBSERVED SPECTRA 

LEs are faint, and therefore they do not contribute 
significantly to the SN spectra around maximum. Fol¬ 
lowing this reasoning, we first analyze the reddening at 
maximum light, to determine the total extinction or Ay 
due t o C SM and/or ISM using standard extinction laws 
(see ^ 4.3). Then, analyzing the best fits to the spectra at 


later epochs we can discern if it is compatible with the 
extinction found at maximum light or if it is necessary to 
include LEs or, alternatively, reduce the amount of dust 
in the line of sight. 


4.1. Bata 

The spectra we use we re taken from the C arnegie 


Supernova Program CSP (Folatelli et al. 
data o f the Center for Astrophysics CtA 


20131), public 


Ilondm et al. 


2012b), the Berkeley Supernova la program BSJNIP (Sil- 


verman et al. 2013b) and The Onhue Supernova Spec- 
trum Archive SUSPECT (see tablej^ . We analyzed only 
the subset of SNe that were already classified as fast or 
slow Lira decliners during the Lira Law phase in F13. 


First, the spectra are co rrected for Milky Way e xtinc- 
tion using the values from Schlafly &: Finkbeiner (2011) 
and deredshifted to rest-frame . Then we smooth the 
spectra using a non-parametric fit with a velocity win¬ 
dow of 1000 km/s and a wavelength regridding of 5 A. 
We also compute the dispersion of each original spec¬ 
trum with respect to its smoothed version to estimate 
the noise in our smoothed spectra. We normalize each 
smoothed spectrum to the same V band flux by numeri¬ 
cally con volving the sp ectrum with the filter transmission 
function (Besselll 19901, in order to put all spectra in the 


same scale and be able to compare the shape and fea¬ 
tures of the spectra instead of the absolute fluxes, which 
are difficult to calibrate precisely. We also adjust the 
shape of the spectra to mat ch the observed co lors inter¬ 
polated to the given epoch (jHsiao et al.||2007|), in order 
to have spectra consistent with the available photometry 
and with the previous work in F13. 

We analyze the different SN spectra at maximum light 
and during the Lira law phase at 5 different epochs or 
time windows centered at 40, 50, 60, 70 and 80 days af¬ 
ter maximum light with a width of 10 days. To have a 
single representative spectrum at every time window per 
SN, we make weighted average spectra with the avail¬ 
able spectra of each SNe. For more details about this 
weighted averages check 


4.2. Template spectra 

F13 showed that slow Lira decliners present weaker 
equivalents width (EW) of blended Na I D1 & D2 narrow 
absorption lines, while fast Lira decliners have stronger 
EWs and redder colors independent of environmental 
factors. One possible interpretation of these results is 
that fast Lira decliners have CSM that produces LEs or 
that they have nearby dust that is sublimated in time. 
In order to test these hypothesis and analyze the spectra 
of fast Lira decliners with our models, we compare them 
with some standard unreddened and non-evolving CSM 
spectral time series given in this case by the slow Lira 


decliners. For this, we construct different template 
spectra that cove r the intrinsic S N la variability at 
different epochs. Chotard et al. (2011) showed that 
most of the intrinsic variability of spectra in SNe la at 
2.5 days after maximum light can be characterized with 
the equivalent width (EW) of Si II 4131 A and Ca II 
H&K. We follow the approach and use these two lines 
and a stretch parametrization to describe the shape of 
the light curve and defined as a factor multiplying the 
time axis ( Perlmutter et al.|1997 Goldhaber et al.|2001 1 
to construct different templates. However, we can not 
measure the EWs of the Ca and Si lines in all of our 
SNe la. Given than the EW of Si II 4131 is correlated 
with the light curve stretch parameter, we decide to use 
epoch and stretch as our main variables to construct 
templates accounting for the intrinsic variability in SN 
la spectra. 


To ensure that these average templates are not heavily 
biased by few extreme SNe, we perform a bootstrap sim¬ 
ulation, i.e., we constructed 100 different templates using 
random sets from the original. For more details see 
In figure]^ we compare the bootstrap templates for fast 
and slow Lira decliner spectra at maximum light and 50 
days later, and one standard deviation (a) regions. At 
maximum light there is a clear difference between fast 
and slow Lira decliners near 4000 A, the average differ¬ 
ence is about 1.4(7 between 3400 - 10000 A and O.Scr if 
we do not consider SN 2006X and SN 2003cg. Between 
3400 - 5500 A the difference is about 2.3a and 1.4cr if 
we do not consider SN 2006X and SN 2003 eg. But this 
difference after 50 days decreases to O.Scr between 3400 


- 10000 A (0.9cr without SN 2006X and SN 2003cg) and 
0.5cr between 3400 - 5500 A (0.6cr without SN2006X and 
SN2003cg). We also find that the dispersion among fast 
Lira decliners is larger than in our slow Lira decliners 
sample. These results confirm what was found in F13, 
and we stress the fact that they are valid irrespective 
of warping the spectra to the observed colors. Finally, 
we point out that the template of th e slow Lira declin- 
ers is quite similar to the template by|Hsiao et al.| (|2007| 
suggesting that this is the group of more “normal” unred- 
dened SNe la. 


4.3. Extinction law at maximum light 

To find an extinction law for each fast Lira decliner at 
maximum light, we use as unreddened reference a slow 
Lira decliner template representing the sa me epoch and 
stretch , and t he extinction law described in Fitzpatrick <^| 
Massa (1990), adopting the mea n values for the param¬ 
eters found in Fitzpatrick (1999) and leaving the visual 
extinction Ay as free parameter. 

We fix the total-to-selective extinction ratio Ry at the 
standard Milky way value of 3.1. A discussion on this 
is presented in ^ For each fast Lira decliner average 
spectrum we fit Ay and a normalization parameter Ci, 
that corrects the fact that our spectra are normalized to 
their V-band flux, minimizing a chi-square function (for 
more details seeO. 

To analyze gloral differences be tween both s ubset s of 
SNe we fit the extinction law of Fitzpatrick (1999) to 
the bootstrap templates of fast Lira decliners using as 
reference the bootstrap templates of slow Lira decliners 
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TABLE 2 

Nearby SNe Ia spectra used in this analysis (besides the data from CSP, CfA and BSNIP) from The Online Supernova 

Spectrum Archive (SUSPECT). 


Name 

SN 2005hk 
SN 1999ac 
SN 1998aq 
SN 2005cf 
SN 2003du 
SN 2005am 
SN 2006X 
SN 1999aa 
SN 2002bo 
SN 2000CX 
SN 1994D 
SN 2003cg 


Sources 


Phillips et al.|2007 Blondin et al.|2012a| Chornock et al.|2006 


Garavini et al.|2005 


Branch et al. 12003 


la ^ 

Garavin i et al.l2007a Wang et al. 12001? ^ Bufano et al. 12009 lLeona rdl2007 
Stanisn ev et ^ |2007|Anupama et al.|2005|Gerardy|2odll * ' 


Stanisnev et ^|2007|Anupama et al.|2005|Gerardy|2005' 


|Leonard|2007 


Yamanaka et al.|2009 _ Wang et aI.|2008b|Sternberg et al.|2011a 


Garavini et al.|2004 


Benetti et al.|2004 

Li et al.|200l' 


||Patat et al. | l9^ _ Gomez L6p ez|l998 

Elias-Rosa et al.|2006| 



Wavelength 'A' 



Fig. 6.— Bootstrap templates for slow and fast Lira decliners 
SNe at maximum light (upper panel) and 50.0 days after maximum 
light (lower panel), for a stretch of 0.98. In blue is the B — V 
slow Lira decliner template , in red (green) is the B — V fast Lira 
decliner template including (without*) SN 2003 cg and SN 2006X. 
In black line a template by |Hsiao et al.| l|2007|l is presented. The 
filled contours represent the one standard deviation regions based 
on the dispersion in the spectra. 

presented in Figure |12[ We find at maximum light an 
extinction with Ay = 0.64 ± 0.04 and Ry = 3.0 ± 0.2 
when we include the highly reddened SN 2006X and SN 
2003cg. When we exclude them we obtain Ay = 0.24 ± 
0.01 and Ry = 3.3 ± 0.1. However, during the Lira law 
we obtain an extinction with Ay = 0.17±0.01 and Ry = 
3.3 ± 0.1 significantly lower than at maximum light and 
when we exclude these SNe we obtain Ay = 0.04 ± 0.01, 
Ry = 3.6 ±0.1. These results recover some of the results 
in F13, fast Lira decliners show a greater extinction at 
maximum light, which starts to decrease during the Lira 
law. Moreover, their Ry slightly increases over time. 

4.4. Light Echo fit 

Using the model described in 0 plus the values for 
the dust albedo w and the degree of forward scattering 
g, corresponding to dust with Ry = 3.1, we fit the LE 
model to observed spectra of fast Lira decliners in the 
Lira phase. The parameters to fit are /csm and a nor¬ 
malization constant C 2 to correct for the fact that our 
spectra were normalized by their flux in the V band. 

The function to minimize for each average spectrum i 


2 ^ V (r(^.A)-C2±ft,/cSM,A))^ 

A* 2^ sf{t,X)^ + {C2SF{t,fcsM,W 

where 

F{t, A) = f{t, A)10"° ‘‘^^^’ -b S{t, (8) 

To evaluate S{t,X) we use eq. |^in a slightly different 
version because we do not have the intrinsic fluxes in 
our spectra as they were previously normalized (for more 
details, see 0. 


5. RESULTS 


5.1. Extinctions 

Out of 31 individual spectra from different fast Lira 
decliners at maximum light, minimizing the of eq. Cl 
we obtain 24 SNe with positive Ay values and 7 with 
nonphysical negatives values. Those SNe are excluded 
from our sample for the LE analysis as they present less 
extinction than the template at maximum light, making 
the posterior fit of our LE model impossible. This may be 
caused by some SNe in the sample of slow Lira decliners 
having non negligible host extinction. At maximum light 
the mean Ay is 0.44 while excluding the highly reddened 
SNe, SN 2006X and SN 2003cg, it decreases to 0.27. But 
during the Lira law the amount of extinction starts to 
decrease to values near zero. In Table |3] the mean dif¬ 
ferences in Ay per SN are presented with respect to the 
value found at maximum light excluding those SNe with 
negative Ay at that epoch. This result confirms that 
fast Lira decliners are more extincted at maximum light 
than slow Lira decliners, and during the Lira law both 
groups become more similar, as shown in F13. This is 
also valid without the warping of the spectra to the ob¬ 
served colors. However, there are a few SNe that show a 
positive A Ay. This could be produced by an artifact in 
our templates or a bad extinction law fit. 

We find a small correlation between the Ay values 
found at maximum light and the B — V slope of the 
Lira law, with a Pearson’s correlation coefficient of -0.5 
(see figure^, which is consistent with the analysis done 
in E13. 
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TABLE 3 

Mean and uncertainty of AAy = Av(t) — Ay (max) for fast Lira decliner SNe at different epochs. 


40 days 

50 days 

60 days 

70 days 

80 days 

-0.26 ±0.07 

-0.19 ±0.09 

-0.54 ± 0.11 

-0.21 ± 0.17 

-0.29 ±0.18 


TABLE 4 

Comparison of ISM, LE and DS models with overall spectra: fraction of fast Lira decliner SNe with positive Av at 
maximum light that have a BIC value which favors the LE or DS models vs the other models at different time windows. 



40 days 

50 days 

60 days 

70 days 

80 days 

Number of SNe 

16 

11 

8 

10 

7 

LE vs ISM 

0.50 

0.36 

0.63 

0.5 

0.57 

LE vs ISM & DS 

0.19 

0.18 

0.13 

0.2 

0.14 

DS vs ISM & Echo 

0.44 

0.27 

0.63 

0.3 

0.43 


TABLE 5 

Comparison of ISM, LE and DS models with specific spectral regions: fraction of fast Lira decliner SNe with positive Ay 
at maximum light that have a BIC value which favors the LE OR DS MODELS vs THE OTHER MODELS AT 50 DAYS AFTER MAXIMUM 

LIGHT FOR DIFFERENT WAVELENGTH REGIONS. 



3900-10000 A 

3000-6500 A 

6500-10000 A 

3600-4800 A 

5900-6400 A 

LE vs ISM 

0.36 

0.45 

0.45 

0.27 

0.18 

LE vs ISM & DS 

0.18 

0.36 

0.45 

0.27 

0.18 

DS vs LE & ISM 

0.27 

0.18 

0.0 

0.0 

0.0 



Fig. 7.— Av fit results at maximum light vs B — V slopes of the 
Lira law from F13. 


5.2. Light echoes 

Once an Ay per fast Lira decliner SN is found, we fit 
our light echo model to the different fast Lira decliner 
spectra to search for spectral evidence of LEs. We look 
for the minimum of eq. Q by varying /csm and C 2 . We 
fix Rv to 3.1 and the ramus of the CSM shell to 0.05 pc 
(more on this in 

We also fit the data with two other models to compare 
the goodness of fit of the light echo model applied to the 
observations. The three models that we compare are: 

1. The light echo model in which we fit for /csm and 
C 2 . It also includes ISM extinction to get a total 


extinction consistent with the amount of extinction 
found at maximum light. 

2. A pure ISM model that just uses a late time spec¬ 
tral template with the same extinction found at 
maximum light. In this scenario we just fit a nor¬ 
malization constant Ci. 

3. A dust sublimation (DS) model consisting of a new 
pure extinction fit during the Lira law, with a lower 
Ay value than the one found at maximum light. 

Then we calculate the Bayesian Information Criterion 
(BIC) parameter in order to compare the goodness of 
the three models applied to the same data. The BIC 
parameter is obtained as 

BIC = + k ■ ln(n) (9) 

where k is the number of free parameters to fit in each 
model and n is the number of data points or wavelengths 
with measured flux values. The best model is the one 
with the lowest BIC. It is a combination of plus a 
function that penalizes having too many free parameters 
overfitting the data. With the BIC parameter we judge 
whether the LE model is able to explain the observed 
spectra better than a simple extinction or DS model. The 
results of the fits to the entire spectra and comparison 
between the three models are shown in Table lU 

According to our results, the LE model works better 
than using just the extinction law found at maximum 
light in 50% of the cases. But when we compare it to the 
dust sublimation model, the fraction of favorable cases 
drops to values lower than 20%, giving larger weight to 
a scenario where Ay simply decreases with time. 

We also fit the three models to the spectra in par¬ 
ticular wavelength ranges: blue region (3000 - 5000 A); 
red region (6500 - 10000 A); signature I (3600 - 4800 A) 
































and signature II (5900 - 6400 A). The results are shown 
in Table The analysis in wavelength ranges is not 
favourable for the DS model, and the simple extinction 
found at maximum light works better than the LE model 
in most cases. However, this can be explained by the fact 
that when we look at particular regions of a spectrum a 
degeneracy appears between the normalization and the 
extinction, as there is not enough wavelength range to 
anchor the reddening, obtaining a better BIC with the 
constant extinction model than with the DS model which 
has more parameters. 

For each late-time SN spectrum, we derive the value for 
the CSM fraction /csm- Multiplying this value with the 
Ay found at maximum light we can infer the fraction of 
the extinction due to CSM under the hypothesis of LE. 
For the four SNe with Ay > I at maximum light, we 
obtain that their Hy’s due to CSM are lower than 0.5, 
validating our single scattering assumption. 

In Figure we show /csm x Ay for different SNe at 
different epo^s. Only epochs in which the LE model was 
the best according to its BIC value are plotted. If the 
CSM were not disturbed by the supernova we would ex¬ 
pect /csm X Ay to stay almost constant for each SN. Un¬ 
fortunately, only one object, SN 2003W, has more than 
one late-time spectrum consistent with LE to allow us to 
perform this test. The three resulting /csm x Ay values 
range between 0.03-0.23, which we consider a satisfac¬ 
tory agreement considering the simplifications of the LE 
model. 
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Fig. 8. — Ay X /csm vs time since maximum light. Only cases 
in which the LE model works better than the pure extinction and 
DS models are shown. The error bars represent three standard 
deviation errors on the parameters that minimize the function 
under our LE model assumption. 

The number of SNe for which the LE model has a favor¬ 
able BIC in comparison with the ISM and DS models is 
small (7 SNe), and just SN2003W appears in more than 
one epoch. Even in such cases, a visual inspection of 
the spectra does not reveal the signatures expected from 
LEs. This poses serious questions on the LE scenario. 

5.3. Line comparison 

Another way to test our LE model is to use the line 
diagnostics presented in section For this we calcu¬ 
late for all our spectra the pW and pseudo-continuum 


slope for the four diagnostic lines previously defined, in 
the same way we measure them in the simulated spectra. 
We obtain these features at all available epochs and also 
normalize the pW curves by their values at maximum 
light in order to study their evolution. Finally we com¬ 
pare the results for the slow and fast Lira decliners with 
the simulated spectra. Apart from some slight differences 
between the two samples, we find that in general the pop¬ 
ulations are consistent within the errors. Both groups of 
SN la show similar trend and dispersion in their evolu¬ 
tion. As an example. Figure]^ shows the distribution of 
pW at 55 days after maximum light for the two SN sam¬ 
ples. In the context of the LE model, we would expect 
both a difference in the dispersion of pWs and their mean 
value, neither of which is observed. 

To measure possible statistical differences we perform 
a K-S test. We find that for line I the two populations, 
slow and fast Lira decliners, have a 99% probability of 
being drawn from the same distribution. We check all of 
the other lines for which we obtain low K-S values, yet we 
do not find any hint of LE signatures among fast Lira de¬ 
cliners according to our predictions (KS values between 
0.2 and 0.99) . The pW distributions are consistent for 
both populations during the Lira phase. However, some 
pseudo-continuum slopes show very low KS values (lower 
than 0.01). This difference in both SN samples can be 
explained by the greater color dispersion of fast Lira de¬ 
cliners which also tend to be redder at maximum light, 
as shown in F13. As a matter of fact, the trend goes 
in the opposite direction of what is predicted by our LE 
models. 



Fig. 9.— Histograms of the ratio of pseudo-equivalent width at 
maximum light and at 55 days past maximum light for the diag¬ 
nostic echo line 1 for the sample of slow (blue) and fast (red) Lira 
decliners. Vertical dashed lines show the median of the population: 
5.4±1.70 for slow and 5.7dil.57 for fast Lira decliners, respectively. 
The KS test for these distributions gives a probability of 99% of 
being drawn from the same distributions 

Although we do not find a correlation between the evo¬ 
lution of lines and fast Lira decliners, our CSM model 
predictions, especially for line I, can be important to 
diagnose the presence of LEs in other samples of SNe la. 

6 . DISCUSSION 
6.1. Extinction laws and Ry 
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To calculate the extinction a t maximum, we also 


use other reddening laws such a s Card elli et al. (19891 
with the incl usion of O’Donn^ (1994), obtaining simi¬ 
lar results toll'itzpatrickl (|l999() when keeping Ry fixed. 
We al so explore the reddening law proposed in |Goobar| 
(2008), but we discard it for two reasons: the reddening 
law becomes degenerate with our normalization constant 
Cl and also because it accounts for the observed redden¬ 
ing in the context of CSM without considering the evo¬ 
lution of the radiated spectrum. In contrast, we aim to 
find an intrinsic reddening law and the time evolution of 
the spectrum caused by LEs. 

To be consistent in our predictions and in our fits we 
used extinction laws with a fixed Ry of 3.1 as we use 
standard MW values for the albedo and phase function 
of interstellar dust. Therefore our analysis is restricted 
to a statistical point of view more than the study of par¬ 
ticular cases. The Ry value can vary depending on the 
properties of the dust (e.g. grain size distribution and 
composition) and seems to be differe nt from the MW in 


the line of sight towards some SNe la (Mandel et al. 2011 


_ sig- _ 

Burns et al.|2014 ). It is not very clear what range of val¬ 


ues are consistenf with circumstellar dust surrounding a 
SN. It is very important to differentiate in extinction laws 
analysis between the intrinsic Ry, which comes from the 
dust properties, and the observed Ry when a pure ex¬ 
tinction law is assumed, omitting more complex interac¬ 
tions. A future improvement in our model is to compute 
and use the specific opacities, albedo and phase function 
given any dust grain size distribution and composition. 


6.2. Light echo models 

We have explored the possibility of detecting LEs due 
to CSM in SN la spectra. Our results show that LEs 
are not a global phenomenon on fast Lira decliner SNe 
during the Lira law phase. Even though we find that for 
^50% of the spectra the LE model works better than the 
extinction law derived at maximum light, the number of 
favorable cases drop to values near 15% when we compare 
them with the DS models. 

We also fit the LE model using a CSM radius R of 
0.01 and 0.25 pc, instead of 0.05pc. With the smallest R 
we obtain a lower fraction of favorable cases for the LE 
model compared to the original results with R — 0.05 pc, 
even lower than 50% when the LE model is compared 
with just the pure extinction model (ISM). On the other 
hand, when we fix i? = 0.25 we recover almost the same 
results than the original CSM scenario. Therefore, if 
CSM is present, larger radii of 0.05-0.25pc are favoured. 
In principle, it is possible to fit at the same time R, fcsM 
and C 2 , but this is computationally expensive and could 
overfit the data. 

Our CSM model consists of an isotropic spherical shell 
in the limit of negligible thickness. We did not consider 
multiple scattering, which at optically thin scenarios is 
negligible. Nevertheless, we know that multiple scatter¬ 
ing could become important at optical wavelengths when 
the optical depth is larger than 1, i.e. Ay > 1. There¬ 
fore we are unable to analyze SNe with an expected Ay 
due to CSM larger than one, but according to our results 
none of the SNe in our sample presented an Ay due to 
CSM larger than 1 (including SN2003cg and SN2006x). 
A model that includes multiple scattering is necessary to 
predict the effect on the light curves and spectra when 


an optically thick CSM is present (|Amanullah & Goobar 
2011|||Patat||2005| ). 


The CSM geometry may probably be different from a 
spherical shell, e.g. non-isotropic disk or ring geometries 
formed from a planetary nebula have been proposed re¬ 


cently to model time variable Na I D absorption (Soker 


2014). The predicted LE signatures might vary depend- 


mg on the CSM geometry and orientation to the ob¬ 


server. 

For an optically thin shell, a rough estimate of the total 
dust mass in our models can be obtained for a given CSM 
radius, an Ay and a typical ISM dust opacity: 


Md = 6.4x10"® 


R 


0.01 pc 



hv 


8.55 X lO^cm^g- 




( 10 ) 


This is a very high mass in form of dust for a CSM. 
However, we are assuming a spherical shell and a spe¬ 
cific opacity corresponding to interstellar extinction in 
the MW. If we consider a different geometry for the CSM 
or a larger specific opacity, i.e. smaller grains, the in¬ 
ferred mass will vary by orders of magnitude. 


6.3. Light echo effects on the light curves 

Usin g Monte Carlo simulations, [Amanullah fc Goobar| 
(2011) found that different radius or E{B — V) for the 
CSiVI could affect the B — V evolution during the Lira law 
phase in different ways. In order to test the hypothesis 
made in F13 that LEs could increase the B — V decline 
rate during the Lira law phase, we investigate whether 
our CSM model affects the behaviour of the light curve 
during the Lira law phase. We find that LEs make the 
opposite, they tend to smooth the color evolution. This 
goes in the opposite direction of our goal of finding light 
echoes in fast Lira decliners spectra and favors the DS 
model. 



Days after maximum 


Fig. 10.— Simulated light curves using our LE model. In black is 
the original light curve without extinction. The reddest curve rep¬ 
resents a light curve with pure extinction (Av^=0.5 and i?y=3.1) 
and the rest are scenarios with different fractions of CSM. 


In Figure 10 we present these simulated light curves. 
As expected, in the pure extinction scenario the B and 
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V light curves (red line) are just uniformly shifted down¬ 
ward with respect to the original light curve without ex¬ 
tinction (black line). On the other hand, the presence 
of LEs due to GSM modify the shape of these curves, 
increasing the brightness in B and V at later epochs 
as blue light from maximum light is reaching the ob¬ 
server. However, the slope of the B — V evolution actu¬ 
ally becomes shallower during the Lira law, see Table 
This can be explained qualitatively as our LE model ados 
light emitted at previous epochs to the intrinsic emitted 
light, including late time emission with small time de¬ 
lays. These contributions make the observed light curve 
evolve slowly. 

The wiggles in the B and B — V light curves are not 
real. They are caused by the way we compute the B 
magnitudes from the template spectra, which depend on 
the available SN spectra at each time window. On the 
other hand, the V magnitudes match the observed pho¬ 
tometry by construction. Despite this, the general shape 
of the light curves is clear. 

TABLE 6 

B — V SLOPE DURING THE LIRA LAW PHASE FROM SIMULATED 
LIGHT CURVES. 


fcSM 

Slope ± error [mag/day] 

0.0 

-0.010 ±0.002 

0.2 

-0.008 ±0.001 

0.4 

-0.008 ±0.001 

0.6 

-0.007 ± 0.001 

0.8 

-0.004 ± 0.001 

1.0 

-0.005 ±0.001 


6.4. Dust sublimation and its effects on the light curve 

A decreasing extinction or opacity could occur if the 
CSM dust that was extincted at maximum light got sub¬ 
limated by an increment of temperature due to the SN 
radiation, reducing the total opacity in the line of sight. 
This sublimation could also change the observed Ry as it 
might change the grain size distribution and composition, 
explaining the evolution found in F13. If this happens 
during the Lira law phase, it will be reflected in a steeper 
B — V slope. To test this hypothesis we simulated li ght 
curves with a decreasing Ay, and variable Ry. Figure ^ 
presents these three different scenarios. We found that 
a decreasing Ay or increasing Ry can make the B — V 
evolution become steeper. 

A smaller Ry makes the extinction law more sensi¬ 
tive to the blue than the red, which is expected if the 
grain size distribution favors small sizes. An increasing 
Ry could be produced if the smaller grains disappear as 
sublimation occurs. Another possibility is that the in¬ 
trinsic CSM Ry was lower than the ISM Ry, therefore 
as the circumstellar dust is being sublimated, the total 
Ry increases reaching values similar to the interstellar 
Ry. This hypothesis is also consistent with the mean 
AAy found at different epochs for fast Lira decliners. 

Wit hin this context , it is important to note that re¬ 
cently Goobar et al. (20141 showed that SN2014J, an 
extincted red SiN la with low Ry dFoley et al.||2014|) 
and w ith strong narrow absorption features (Welty et al. 
20141, has possible signatures of cooling from shocked 


material from nearby CSM of dimensions larger than 1 



-20 0 20 40 60 80 


Days after maximum 


Fig. 11. — Simulated light curves using extinction laws that vary 
over time. In black is the original light curve without extinction. 
The blue curve is a model with Av decreasing from 0.5 to 0.0 and 
Ry = 3.1 held constant. The red curve represents a light curve 
with Av=0.b and Ry evolving from 2.0 to 3.1. The green curve 
represents a light curve with Ay=0.b and Ry evolving from 3.1 to 
2 . 0 . 


Rq. But the scales of this CSM are smaller than the 
ones considered in this paper. 

Another possibility is that the radiation pressure (RP) 
of the SN is blowing away the CSM dust particles. If 
the CSM radius is increased, the observed extinction will 
decrease as the column density decreases, also producing 
blueshifted absorption lines. A rough estimation of the 
time scales in which the RP could produce an observable 
change in the CSM extinction can be calculated as 


trp = 


trp = 


Rcsm 

Av 

md Rcsm c 

fsN CFa At 


( 11 ) 

( 12 ) 


where R is CSM radius, Av is a characteristic velocity of 
the dust grains after they absorb linear momentum from 
the SN incident radiation flux /sn that can be calculated 
using the well known relation between the radiation pres¬ 
sure and flux of an electromagnetic wave. We consider 
spherical dust particles of radius a, internal density of 
I gr/cm^, mass rud, and with a cross section a a calcu¬ 
lated using the Mie theory. Assuming a typical SN la 
luminosity and a time range At= 25 days, centered at 
maximum light, in which the RP injects momentum to 
the dust particles, we can obtain an estimate of the time 
scale in which this effect could be observed: 


Thus, if the radiation pressure is the cause of the de¬ 
creasing extinction at trp ~ 80 days past explosion (Lira 
law), the CSM dust particles should be smaller than 1 fim 
or be at distances of ~ 10^ astronomical units. At these 
distances we expect that the sublimation time scales of 
these smaller grains are much shorter than the calculated 
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above. 

Therefore, even if dust sublimation time scales are too 
long to account for the changes in extinction and Ry, 
during the Lira law phase, we expect that the CSM ex¬ 
pansion due to RP should increase the observed Ry and 
decrease the extinction, as smaller grains are blown away. 

An alternative scenario to the CSM sublimation and 
RP that could explain the decrease in extinction is the 
transversal expansion of the ejecta in a non-homogeneous 
ISM. This possibility is explored in F13 where they con¬ 
clude that this scenario could explain the change in the 
average column density as the photospheric radius in¬ 
creases, but it does not explain the change in the ob¬ 
served Ry over time. 


7. CONCLUSIONS 

A CSM model producing light echoes (LEs) has been 
developed. It is simple enough to be computed quickly 
and be used in our fitting routines, but with enough com¬ 
plexity to account for the albedo of the dust and the scat¬ 
tering phase function. Our model predicts two spectral 
signatures produced by LEs at 4100 and 6200 A that can 
help us discriminate between a pure extinction or extinc- 
tion-|-LEs scenario. These features appear within small 
wavelength ranges as opposed to overall color changes 
that can also be produced by reddening. The evolution 
of these signatures is another tracer of the presence of 
CSM producing LEs. 

We compare our models with observed SN spectra. We 
find that LEs from CSM at 0.01-0.25pc are not a global 
phenomenon in fast Lira decliner SNe when they are com¬ 
pared to slow Lira decliners. ISM or CSM dust being 
sublimated (DS) at later epochs explains better the ob¬ 
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served spectra when the models are fitted using the over¬ 
all spectrum. Additionally, we find no evidence for LEs 
based on the narrow spectral diagnostics predicted by 
our model. 

We explore the effects on the light curve of circum- 
stellar dust being sublimated or blown away by radia¬ 
tion pressure, finding that both scenarios could produce 
a faster B — V decline during the Lira law and a change 
in Ry, although a more rigorous physical modeling is 
needed to explore these possibilities. 

We laid out several ways to improve our models: 
adding a Ry as a free parameter in our fits; using differ¬ 
ent CSM geometries to test our predictions; and a Monte 
Carlo radiative transfer simulation to see if the LE sig¬ 
natures remain in the multiple scattering scenario. The 
analysis can also be improved using a larger sample of 
SNe, particularly those highly extincted and with possi¬ 
ble CSM characteristics. 
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APPENDIX 


AVERAGE SPECTRA PER SN PER TIME RANGE 


For the jth spectrum of the ith SN we defined the following weights: 


w) = 


(Al) 


where to is the epoch that we want to represent and at = 5.0 days is half of the time window. The average spectrum 
for a particular SN at the given epoch to was then calculated according to 


r(A) 




(A2) 


where /j(A) and Sfj{X) are the normalized flux of the jth spectrum and its error. The sum goes over all the available 
spectra of the ith SN with measured fluxes in the given time range. We also computed an error and a representative 
epoch for each average spectrum in a similar way. 


TEMPLATE SPECTRA 


To construct template spectra we define two weight factors for time and stretch of the ith average spectrum of each 
slow Lira decliner SN: 




(Bl) 

(B2) 


We chose at = 1.5 days and CTs = 0.11 in order to reproduce a specific epoch and stretch. We do not have strong 
arguments to choose particular values, thus we use the standard deviation of stretchs in our sample for as and a at 
value smaller than our time windows, but large enough to reproduce smooth light curves. Defining at = wl x wl the 
template spectrum with a certain epoch to and stretch sq is 


EiQ^i/»(A)M/z(A)^ 
E* a,/(5/,(A)2 


(B3) 


where the sum goes over all the available already averaged spectra from slow Lira decliners with measured fluxes at 
wavelength A. Finally we normalized the template spectrum by its flux in the V band. 

To ensure that these average templates are not heavily biased by few extreme SNe, we performed a bootstrap 
simulation. Then computing the mean template and the dispersion around it, we obtained a “bootstrap template” 
and its error. In figure we compare a template constructed using weighted averages and a template using the 
bootstrap simulation. In both cases we used all the slow Lira decliners sample. There are slight differences between 
both spectra, but these are not significant and the estimated errors are very similar. 


EXTINCTION LAW FIT 

To fit an extinction law we minimize a chi-square function for each average spectrum i 

.,2 v (.r(A)-/°(A)I0-°-4^(^)M)^ 

* A + 

where /®(A) is the normalized flux of the ith fast Lir a de cliner SN at maximum light and f^{X) is an unreddened 
template representing the same intrinsic flux. In figure we present a particular extinction law fit for SN1999gd at 
maximum light. We found that Ay = 1.280 ± 0.003 best represents the extinction using our template spectrum. 


LIGHT ECHO FIT 


To fit our light echo model we introduce a new factor N to force the spectra to evolve consistently with the light 
curve in the V band. Hence we calculate S{t, A) as 

S{t) = /o"”“ T)f{t - r)<E>'(r)dr (Dl) 

N{t,T) = (D2) 

where we omitted the wavelength dependency of S, f^, ui, and A (the extinction law found at maximum light). 
f^(t) is a template spectrum representing the same intrinsic flux at time t. We multiply each spectrum by Nit^r), 
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Fig. 12.— Template spectrum of a slow Lira decliner SN using weighted averages vs using the bootstrap technique at 50 days after 
maximum light and using a stretch of 0.98. 



Fig. 13.— Extinction law fit for SN1999gd at maximum light using|Fitzpatrick| (|1999|l extinction law. The best fit with Ry = 3.1 yields 
Ay = 1.280di 0.003. In red we present the observed spectrum of SNlfiyygd at maximum light, in black a constructed template representing 
the same epoch and stretch. The blue line is the dereddened spectrum and in green, the extinction factor or 10“*^'^^^ shape. 

which is the V light curve normalized in to correct for the fact that all our templates are normalized by their flux in 
the V band. This factor is computed using weighted averages of the V magnitudes of slow Lira decliners considering 
epoch and stretch a s we did for the template spectra. The V magnitude data was taken from SiFTO fits to the data 
( Conley et al.||2Q08 ). 

































